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1 Introduction 



The perturbative series for QCD observables is not convergent, the n -order coefficients ex- 
hibiting n\ growth at large n. In the Borel plane there are singularities along the real axis: 
ultraviolet (UV) renormalons along the negative real semi- axis, and infrared (IR) renormalons 
lying on the integration contour, along the positive real semi-axis. The latter render the Borel 
integral ambiguous, the ambiguity being structurally identical to ambiguities arising in the 
non-perturbative operator product expansion (OPE) [lj. Consequently, all-orders perturba- 
tion theory is only well-defined if supplemented by the non-perturbative OPE, allowing the 
ambiguities to cancel. 

In practice all-orders calculations are only possible in the "leading-6 approximation", 
in which the n th -order perturbative coefficient is recast as an expansion in powers of b = 
(33 — 2Nf)/6, the first beta-function coefficient (in SU(3) QCD, with Nf active quark flavours). 
The "leading-6" term (proportional to 6 n ) is then used to approximate the n th -order coefficient 
[21E1H]. Crucially this term can be identified to all-orders from large Nf calculations, involving 
a restricted set of diagrams in which a chain of n fermion bubbles is inserted into a basic 
"skeleton" diagram, in all possible ways. In this way one can perform an all-orders leading-6 
resummation by explicitly evaluating the Borel transform, and defining the Borel sum, suitably 
regulated to control the IR renormalon ambiguities. Such results have been used in the past 
to try to assess the likely accuracy of fixed-order perturbation theory for quantities such as 
the R e + e - ratio at low values of the cm. energy Q 2 (where the QCD coupling is not so small), 
and for its tau-decay analogue R T (where again, at the scale of the r mass, the coupling is 
relatively large) [2j EJ El El [8] . DIS sum rules have also been studied in this context [9l [TUl fTTl 

mi ei ng. 

In matching to the exact fixed-order perturbative calculations (NLO and NNLO coeffi- 
cients are available for R e + e -, R T and DIS sum rules), one faces the problem that the all-orders 
leading-6 result is only renormalization scheme (RS) invariant if the one-loop form of the cou- 
pling is used, whereas the exact result must involve the higher-loop coupling. The matched 
leading-6 resummed result consequently depends on the assumed renormalization scale at which 
the matching takes place. To avoid this ambiguity it has been argued in the past that one 
should remove scale dependence by a resummation of scale logs to all-orders; employing the ef- 
fective charge approach [16], or the "complete renormalization group improvement" (CORGI) 
approach [17] . These approaches involve scheme-invariants which can be approximated at the 
leading-6 level and resummed [71 [18] . 
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In this paper we will use the CORGI approach to perform leading-6 resummations for 
DIS sum rules. The aim will be to assess the reliability of fixed-order NNLO perturbation 
theory for these quantities, at the rather small values of Q 2 for which they have been mea- 
sured [191 [2DJ |2TJ [251 ESl IS] • I n addition we shall used a recently proposed model for power 
corrections, based on the QCD skeleton expansion [25], to obtain fits for the non-perturbative 
contributions. 



2 Fixed-order and all-orders predictions 

The three physical quantities we are interested in are the GLS sum rule [26], the polarized 
Bjorken sum rule [27], and the unpolarized Bjorken sum rule [28]. These quantities can be 
written as a parton model result plus perturbative corrections (denoted by calligraphic letters), 
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(3) 



where |<M/flv| = 1-257 ± 0.003 [29], and in the case of Kgls{Q 2 ) we have neglected contribu- 
tions from 'light-by-light' scattering diagrams. The perturbative corrections take the form of 
an expansion in powers of the strong coupling a = a s /ir, 



K(a) =a(Q 2 ) + J2k n a n+1 (Q 2 ), 



U{a)=a{Q 2 ) + Y J u n a n+ \Q 2 ). (4) 

n=l n=l 

These expressions can be evaluated at fixed-order, and the coefficients k n and u n are known 
up to n = 2 [301 E] (NNLO). 



^nnlo(o) = a + kia 2 + k 2 o? 



^nnlo(o) = a + uia 2 + u 2 a 3 . 



(5) 
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The leading-6 approximation 

The leading-6 approximation arises from rearranging the expansion as one in powers of Nf. 
We can write the n th -order coefficients in the form, 

k n = £>] Nf + k^ 1] N]- 1 + . . . + A4°] , (6) 

In] 

and similarly for u n . The leading Nf coefficient kn can be computed exactly to all-orders as 
noted above. One can replace Nf by (33/2 — 36) to obtain an expansion in powers of b, 

k n — k^ ^ b + k^ 

The leading-6 term k n L ^ = k n n ^b n can then be used to approximate k n [2]. To see how well 
this approximation works we can compare the known exact results for k\, &2, and u±, 112, with 
the leading-6 approximations k[ L \ k^ and u[ L \ 112 ■ We shall assume the MS scheme with 
renormalization scale /j, 2 = Q 2 . For /C(a) one has the exact coefficients 

k x = -0.333iV7 + 4.58 (8) 
k 2 = 0.177^-7.61^/ +41.4. (9) 

To be compared with the leading-6 approximations 

k[ L) = -0.333iV7 + 5.5 (10) 
4 L) = 0.1777V) - 5.86iV7 + 48.3. (11) 

The leading- Nf coefficients of course agree exactly by construction, but one sees agreement in 
sign and ~ 20 — 30% level agreement in magnitude for the sub-leading coefficients as well. The 
corresponding coefficients for U{a) are 

Ul = -0.444^ + 5.75 (12) 
u 2 = 0.2397V) - 9.5iV/ + 54.2. (13) 

To be compared with the leading-6 approximations 

U [ L) = -0.444iV/ + 7.33 (14) 



u 



(L) 
2 



0.239iV) - 7.897V/ + 65.1. (15) 



Again one sees that sub-leading coefficients are well-reproduced in sign and magnitude. This 
can be partially understood since it is possible to derive a fully normalized asymptotic result 
for the k^- n ~ r ^ and iJ n_r ] coefficients for large n and fixed r. Such a result for the QCD Adler 
function of vacuum polarization, D(Q 2 ), was derived in Ref. [32] from a consideration of 
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the four-fermion operators controlling ultraviolet renormalon singularities. It can be easily 
extended to the DIS sum rules since the same dominant four-fermion operator is involved |33j. 
If one writes the leading-6 term as 

*£ L > = ^; + fc%r 1 + ... + ^ 0] (16) 

uSP = u^N] + u^N^ + ... + ut\ (17) 
then one has the asymptotic result for large n and fixed r 

k [ r r] « k [ r r] (i+o(l))+o(±) as) 



n 



u 



[ n _ r] ^ ^-,]( 1 + (I^ +0 ^_1). (19) 



The 0(1/ N 2 ) are terms sub-leading in the number of colours, TV. The asymptotics agree up to 
0(1/ n) terms in the "planar approximation" [32], where at each order in Nf only the highest 
power of N is retained. Of course this asymptotic result does not explain the ~ 20 — 30% level 
agreement observed on comparing Eqs. (8), (9), (12) and (13) with Eqs. (10), (11), (14) and 
(15), for such small values of n = 1, 2. The agreement is also remarkably good for the kn and 
u"n coefficients which correspond to Nf = 0, or the large- N limit. Similar remarks apply for 
the Adler function D(Q 2 ) and some speculative ideas as to how this might arise are given in 
Ref. [32J. Whilst this good agreement for given powers in the Nf expansion is interesting, of 
more relevance is how well the overall k\, hi and u%, v>2 coefficients are approximated, since 
one would not expect the ~ 20 — 30% accuracy to survive the addition of the terms. We give 
below the exact and leading-6 coefficients for Nf = {0, 1, 2, 3, 4, 5}. For K(a) one has 

ki = {4.58,4.25,3.92,3.58,3.25,2.92} (20) 

k[ L) = {5.5,5.17,4.83,4.5,4.17,3.83} (21) 

k 2 = {41.4,34,26.9,20.2,13.9,7.84}. (22) 

k { 2 L) = {48.3,42.6,37.3,32.3,27.7,23.5}. (23) 



Whilst for U(a) one has 



ui = {5.75,5.31,4.86,4.42,3.97,3.53} (24) 



{7.33, 6.89, 6.44, 6, 5.56, 5.11} (25) 



u 2 = {54.2,45,36.2,27.9,20.1,12.7} (26) 



{65.1, 57.5, 50.3, 43.6, 37.4, 31.6}. (27) 



As expected from the above discussion one finds good ~ 20% level agreement for Nf = which 
gradually worsens to factor of 2 level agreement at Nf = 5. 
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All-orders predictions 



The resummation of leading-6 terms to all-orders, which we shall denote as /Coo (a), can be 
accomplished using the Borel sum method, 



n=0 



nl 



/Coo (a) ^ / dze- z l a B\KS L \z). (28) 
Jo 



The Borel transforms of the three above quantities in the leading-6 approximation can be 
calculated from the results in Refs. [3l] and [35]. One finds, 
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1 
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(29) 



(30) 



22 y v z iJ V z 2 , 
The Borel transforms have simple pole singularities on the positive z-axis (IR renormalons), 
and on the negative z-axis (UV renormalons); the positions are at z = z n = 2n/b. One can 
easily compute the PV regulated leading-6 Borel sums, /Coo (a) and U^a) [2], 

1 



/Coo (a) 



9b 



8e 2i/a(Q 2 ) Ei (_ Zl/a (g2 )) + 2e 2 2MQ 2 ) E i(-z 2 /a(Q 2 )) 



+ 16e~ 2l/a(Q2) Ei(z 1 /a(Q 2 )) - We~ Z2/a{Q2) Ei(z 2 /a(Q 2 )) 



(31) 



and, 



Uoo{a) 



1 

36 



-*i/a(Q ) E i(z 1 / a (Q 2 )) _6e^ 2/a(Q } Ei(z 2 /a(Q 2 )) 
2e" 2/a(Q2) Ei(-z 2 /a(Q 2 )) . 



(32) 



(33) 



Ei(x) is the exponential integral function defined (for x < 0, UV renormalons) as, 

/DO g — * 

dt — , 
-x t 

and for x > 0, by taking the PV of the integral (IR renormalons). For these all-orders resum- 
mations to be RS (scale)-invariant we need to take a(Q 2 ) as the one-loop coupling, 



a(Q 2 



61n(Q 2 /A 2 )' 



(34) 



The all-orders results of Eqs. (|3ip and (|32p are derived in the so-called U-scheme 
(MS with renormalization scale n 2 = e _5//3 Q 2 ), and hence A in Eq. (|34j) will refer to that 
in the U-scheme defined by, 

A v = e 5 / 6 A^. (35) 
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CORGI perturbation theory 

As we noted above there is a problem if we wish to match these leading-6 resummations with 
exact higher-order NLO and NNLO calculations for IC(a) and U(a), which involve the higher- 
loop coupling. We shall avoid the matching ambiguity by employing the CORGI approach 
[T7] . The standard RS-dependent coupling a(/i) satisfies the RG-equation, 

= (3(a) = -ba 2 {\ + ca + c 2 a 2 + c 3 a 3 + •••)• (36) 

Here b = (33 — 2Aj)/6 and c = (153— 19Nf)/12b are universal and RS-invariant, the higher co- 
efficients, c n >2, are RS-dependent and may be used, together with dimensional transmutation 
parameter, A, to label the scheme. The CORGI approach consists of resumming to all-orders 
the RG-predictable terms available to a given fixed-order of calculation. The standard pertur- 
bative expansion is then replaced with the following expression, 

oo 

K(ao) = a + ]T X n c% + \ (37) 

n=2 

which removes the renormalization scale dependence of Eq. (|4|) . We will use a tilde to denote 
CORGI-ized results. Initially we shall deal with IC(a) but the results are easily generalized to 
U{a). 

The coefficients X n are renormalization scheme invariant quantities, each of which can 
be derived from an N n LO calculation of the coefficients of /C(a), and of the beta function 
equation coefficients, c n . For n = 2 and 3 they have the form, 

X 2 = k 2 - k\ - ck x + Q2, (38) 

c/c 2 1 

X 3 = k 3 - 3kik 2 + 2kf + -± - hc 2 + -c 3 . (39) 
The coupling ao used in Eq. ([371) has the form [361 EZ] ; 

«»(<?) = ~\ /0 *\-v^V {40) 

and is the solution to the beta function equation in the 't Hooft scheme [38] (a scheme where 
c n >2 = 0). Here W is the Lambert W function, defined implicitly by W (z)exp(W (z)) = z. The 
"—1" subscript on W denotes the branch of the Lambert W function required for asymptotic 
freedom, the nomenclature being that of Ref. |39j . 

As a result of adopting the CORGI approach, we use a new scale parameter A^; which is 
RS invariant but dependent on the observable under consideration. It is related to the MS scale 
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parameter by, 



Ak = A^exp^/b) 

= Am(j) ^eMkr/b). (41) 

Here k^ s denotes the NLO perturbative coefficient with scale choice ^ = Q. The exponential 
term in the above equation has the effect of resumming a set of RG predictable terms present 
in the full perturbative expansion, it also reproduces the 0(a 2 ) term present in Eq. (J3|) but 
absent from Eq. (|37p . The difference between A and A is due to different conventions for 
integrating the beta-function equation, A being the standard convention, and A being the 
convention favoured in Ref. |40|. 



Equation ([37j) can easily be adapted to provide a prediction for the unpolarized Bjorken 
sum rule of Eq. ([3]). We simply use the equivalent coefficients of U(a) in Eqs. (|38l) and (f39l) 
and the following coupling, 

MQ) = ~, f ~\ t e>\-V*\\ ' (42) 

where, 

A u = A m exp(</b). (43) 



Equation (131]) can be modified so that is resums the leading-6 components of Eq. (|37j) 
and essentially becomes an all-orders CORGI result. We define the coupling a v through, 

1 = 1 + ^ 6, (44) 
a v do 

where is calculated in the V-scheme. The coupling ao in Eq. (|44p is that of Eq. (|40p . The 
all-orders CORGI result can now be obtained from Eq. ([3T|) using Eq. 



n=2 

= /Coo(a„), (45) 

where Xn , in analogy with Eq. ([7]), is the leading-6 part of X n , for example: 

Af> = k?-(k?)\ (46) 
X? = kf ) -Sk?k?+2(k?) 3 . (47) 
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We can also improve upon ACoo(cio) by adding to Eq. (|4"5]) the known sub-leading-& part of X 2 , 

^oo+K) = £oc(a ) + (X 2 - xf ] b 2 )al{Q). (48) 

Equation (|48p now contains all the information we have about the perturbative coefficients for 
/C(a) at all-orders. Analogous expressions also hold for Woo (a) and U 00+ {a). 



In Ref. [25], we noted that the one chain result of Eq. (|3ip . is finite at the Landau pole 
(Q = A), and remains so for values of Q below A. Remarkably, this finiteness at Q = A holds 
when we use the 't Hooft coupling of Eq. (|40p . Furthermore, both of these results have the 
same values at their respective Landau poles, i.e. 

K(a) =/C(a ) = -4 In 2. (49) 

v ; <2=A v u; Q=h K 9b v ; 

Similar relations apply to U{a) and also to the Adler-D function. This conclusion is altered 
slightly for the actual CORGI result of Eq. ([43]) because of the use of a v {Q). In this case the 
result remains finite at Q = A^, but has a different value to that in Eq. (|49p 



We can also consider NLO and NNLO CORGI, fixed-order approximations which can 
then be compared with all-orders resummations, 

/C NLO (ao) = ao(Q), (50) 
^NNLo(ao) = a (Q) + X 2 al(Q). (51) 

By comparing these with the all-orders predictions, we can asses the reliability of fixed-order 
perturbative predictions. Before doing so we consider the accuracy of the leading-6 approxi- 
mation for the CORGI invariants X 2 ()C) and X2(U), as we did for the perturbative coefficients 
earlier. For X 2 (/C) one finds that 



bX 2 (K) = -0.0221^1 + 1.54^-29.1^/ + 98.6, (52) 
if + 1.09A^ 



bx[ L \lC) = -0.0221Ar| + 1.09A^-18.1iV7 + 99.4. (53) 



Notice that the invariants expanded in powers of b contain a ft -1 term |41| and so it is necessary 
to multiply by a factor of b to obtain the expansion in powers of Nf. As we found for the 
perturbative coefficients in Eqs. (20-27) there is reasonably good agreement in sign and magni- 
tude of the coefficients with astonishingly good ~ 1% level agreement for the Nf = large- N 
coefficient in this case. The numerical values of the A 2 (/C) invariants for Nf{0, 1, 2, 3, 4, 5} are 

X 2 (K) = {17.9,13.7,9.61,5.48,1.33,-2.87}, (54) 
X ( 2 L) {1C) = {18.1,15.9,14,12.1,10.4,8.78}. (55) 
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We see good agreement for Nf = which gradually worsens for increasing Nf. Note that 
X2(tC) changes sign between Nf = 4 and Nf = 5, and so big cancellations are involved in this 
region. For X 2 (JA) we find 

bX 2 {U) = -0.0139JVJ + 1.18JVJ - 25.1JV) + 87.8, (56) 
bX ( 2 L) {U) = -0.0139JVJ + 0.688iV| - 11.3JV/ + 62.4. (57) 

The agreement appears to be significantly worse in this case although the pattern of signs is 
reproduced. The numerical values are 

X 2 (U) = {16,12.3,8.7,5.01,1.22,-2.69}, (58) 
X { 2 L) {U) = {11.3,10,8.76,7.59,6.51,5.51}. (59) 



3 Comparison of fixed-order and all-orders predictions 

The data available for Kgls(Q 2 ) an d K u Bj{Q 2 ) span a range of energy which includes the 
bottom quark mass threshold. It is therefore necessary for us to choose a method for evolving 
the above expressions through this threshold. We adopt the approach detailed in [42J. At 
a particular energy scale we treat all quarks with masses less than that scale as 'active' but 
massless and we ignore quarks with masses greater than that scale. As a consequence, the 
coefficients k n and u n are now iV/-dependent and hence they will depend on the cm. energy 
scale, Q. 

We also perform matching of the coupling at Q 2 = m 2 . At LO and NLO this amounts 
to demanding continuity of the coupling at the threshold but at NNLO and beyond, this 
continuity is violated. This matching forces us to adopt different values of the scale parameter 
in different Nf regions. This is governed by the following equations, 

h N f 

2 1 

A 2 A 2 f m Nf + l \ bNf+1 ^ /*NLoH(NLo\ 

Vn = A ^V"A^~J X6XP l 2b»f* J' (60) 

where S NLO and 5 NNLO can be obtained from the results in Ref. [42]. Here, Ajv^ is the scale 
parameter in the region where Nf quarks are active, m^, is the pole mass of the / quark and 
b Nf is simply b evaluated for Nf quark flavours. 

We must be careful how we apply this matching to the different results we have obtained. 
Equation (|60p is an MS result and hence we carry out the matching for A^jg and then convert to 
the various other scales we have defined in Eqs. ([35]) . (|4"T]) and (|4*3j) . /C NLO (o) is an NLO result 
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Figure 1: Various predictions for the GLS sum rule, superimposed on to experimental data. 

and hence, when applying Eq. ([60]) to this result, we omit <5 NNLO - The results in Eqs. ([4"5]) . ((15]) 
and (j5"T]) are all, at least, NNLO and hence we use the full result from Eq. ([60]) . 

We wish to compare these results with relevant experimental data. Data is available for 
the GLS sum rule [19] and the polarized Bjorken sum rule j20|. I21| , 22, 23J, where the points we 
plot are those arising from the analysis of [23]. We plot the data in Figs. Q] and [2] along with 
our predictions for the full observables. No experiment has so far measured the unpolarized 
Bjorken sum rule, however the possibility exists that it may be extracted from experiments at 
a future neutrino factory [33]. We have taken = 207 MeV, corresponding to the world 
average value a s (M z ) = 0.1176 [25] . 

Within the considerable error bars we see that the different versions of CORGI all fit 
the data tolerably well. The important conclusion is that below Q 2 ~ 5 GeV 2 the different 
approaches give drastically different predictions. This indicates that fixed-order perturbation 
theory cannot be trusted for these lower energies. Even though we have no experimental data 
points for U U Bj(Q 2 ) we have plotted the different versions of CORGI in Fig. [3J We see that 
fixed-order perturbation theory cannot be trusted below Q 2 ~ 2 GeV 2 . A surprising feature 
evident from Figs. 1 - 3 is that even at Q 2 = 2GeV 2 the fixed-order NLO and NNLO results 
remain extremely close, leading to the superficial conclusion that fixed-order perturbation 
theory is working well. In fact as can be seen from Eqs. (54,56) the CORGI invariant X2 for 
both /C and U changes sign between Nf = 4 and Nf = 5 and so fortuitously is very small 
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Figure 2: Various predictions for the polarized Bjorken sum rule, superimposed on to experi- 
mental data. 




Figure 3: Various predictions for the unpolarized Bjorken sum rule. 



12 



K G ls(Q 2 ) 




Figure 4: All-orders leading-6 prediction for the GLS sum rule. Here we use the MS result of 
Eq. (|64p with renormalization scales of // = 2Q, Q and Q/2. 

in this region. When Xg and higher invariants are resummed the potential unreliability of 
fixed-order perturbation theory is revealed. 

We can contrast these CORGI resummations with a leading-6 resummation in the 
MS scheme. If we choose the scale \i = xQ we have, analogous to Eq. (|4"5|) . 

oo 

K%*(x) = a(xQ) + J2kt ] (x)b n a n+1 (xQ) (61) 

?1=1 

= /Coo^x)) . (62) 

Here a(xQ) denotes the full higher-loop (three-loop if matching to NNLO) MS coupling at scale 
fi = xQ, and the k n (x) coefficients are the MS coefficients with scale /i = xQ. The coupling 
a v (x) is then defined by, 

-4-t = —F - b ( lnx+ - ) . (63) 
a v (x) a{xQ) \ 6J K J 

We can match this all-orders resummation to the exact NLO and NNLO perturbative coeffi- 
cients, obtaining, in analogy with Eq. (|48p . 

K™ + (x) = K™(x) + (kt(x) - bkP(x))a 2 (xQ) + (k 2 (x) - b 2 kf' \x))a 3 (xQ), (64) 

where in this a(Q) is the approximated three-loop coupling, given in [29|. 

The resulting prediction is plotted in Fig. H] for the GLS sum rule. Three different 
matching scales corresponding to x = 2, 1, ^ were chosen. As can be seen, the matched re- 



13 



summed perturbative result is hopelessly x-dependent. The CORGI result IC 00+ (a) corresponds 
to the choice x = e~ kl ^^ b . However, we note that no NLO matching is required in the CORGI 
approach, since an infinite set of RG-predictable terms involving k\ have been resummed to 
all-orders, yielding a ^-independent result. 



4 Non-perturbative corrections 

In addition to perturbative corrections it is expected that there will be non-perturbative, 
Higher-Twist (HT) corrections to the sum rules. In a recent paper [25] we have studied the 
infra-red freezing of Euclidean observables, including the sum rules, in the language of the 
one-chain term of the skeleton expansion. The leading-6 resummations we have discussed can 
be written in the form, 

pi poo 

JCoo(Q 2 ) = j Q dt uj^{t)a{tQ 2 ) + J dtuol v {t)a(tQ 2 ) . (65) 

Here ujjc(t) is the characteristic function which is piecewise continuous at t = 1. The two terms 
respectively reproduce the IR renormalon and UV renormalon contributions in the Borel sum, 
but the skeleton expansion result is also defined for Q 2 < A 2 where the standard Borel repre- 
sentation breaks down. Remarkably, there is continuity and finiteness at Q 2 = A 2 . By consid- 
ering the compensation of ambiguities between perturbative and non-perturbative corrections 
alluded to above we are led to an expression for HT corrections in terms of the characteristic 
function, 

A 2 /A 2 \ 

£ht(Q 2 ) = f . (66) 

Here k is a single overall non-perturbative constant. This result is obtained from the V-scheme 
result of Eq. (|3ip and hence the appropriate scale parameter for Eq. ()66[) is Ay. A similar 
expression holds for U HT (Q 2 ) in terms of ivjf(t). 

The characteristic functions cj;c(t) and LOu(t) can be determined from their Borel trans- 
forms, using a result from [25]. wye(t) has the form, 

^"(t) = §-§*, (67) 



and the equivalent expression for U is, 



u i m = |-*, (69) 

<(t) = ^3- (70) 
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X pBj 

d.o.f. 




-0.166 


0.1007 ±0.12 


1.256/4 


-0.06255 


0.002647 ± 0.011 


1.173/11 


^-NNLO 


-0.217 


0.1313 ±0.12 


1.238/4 


-0.1085 


0.004593 ± 0.011 


1.162/11 




-0.33 


0.1997 ±0.12 


1.804/4 


-0.2823 


0.01195 ±0.011 


1.456/11 


K-OQ + 


0.0216 


-0.01311 ±0.12 


2.089/4 


0.03957 


-0.001674 ±0.011 


1.662/11 



Table 1: Fitted values of the non-perturbative constants kqls an d n p Bj, together with their 
respective x 2 P er degree of freedom. 




Figure 5: Predictions for the GLS sum rule, including both perturbative and non-perturbative 
corrections, fitted to the data by varying the parameter kqls- 



In Figs. [5] and [6] we take the four expressions in Eqs. (f45l) . (f48l) . (|50"1) and (f5T1) . supple- 
mented by the non-perturbative term in Eq. (|66p . and perform fitting to experimental data 

(5) 

for both the GLS and the polarized Bjorken sum rules separately. We fix = 207 MeV, as 
before, and use x 2 fitting to obtain the optimal value of the non-perturbative parameters in 
each case, kqls an d n p Bj- The fitted parameters are summarized in Table 1. We note that 
the coefficient of a 1/Q 2 (twist-4) power correction to either Eq. dU) or ([2]) corresponds to a 
value of — 8AyKGLs/3 or — ^\gA/dv l^-v^pBj, respectively. The corresponding values in GeV 2 
are presented in Table 1. For comparison the central values resulting from a three-point func- 
tion QCD sum rules fit [1H US] are —0.294 and —0.013, respectively. We see that the power 
corrections in Table 1 resulting from fitting to the resummed /Coo+ (a) all-orders resummations 
are significantly smaller, although the errors are large. Connections between power corrections 
for the three DIS sum rules have also been explored in |46j . 
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Figure 6: Predictions for the polarized Bjorken sum rule, including both perturbative and 
non-perturbative corrections, fitted to the data by varying the parameter n p Bj- 



5 Conclusions 

By comparing leading-6 resummations based on the CORGI approach we attempted to infer the 
validity of fixed-order perturbation theory for the GLS and pBj DIS sum rules at small energy 
scales Q 2 ~ 3 GeV 2 , where data has previously been used to extract a s (Mz) [29] ■ Figures 1 
and 2 indicate that fixed-order perturbation theory and resummed all-orders predictions start 
to differ drastically for energies below Q 2 ~ 5 GeV 2 . The use of the CORGI approach ensures 
that all RG-predictable scale-dependent logarithms are resummed to all-orders, thus avoiding 
the need for NLO matching which would otherwise make the conclusions about the validity 
of fixed-order perturbation theory dependent on the chosen matching scale (as illustrated in 
Fig. 4). 

We performed fits to data for power corrections using the model proposed in Ref. [25] 
based on the skeleton expansion characteristic function. The size of power corrections in- 
ferred by fitting the all-orders CORGI resummations was much smaller than the central values 
obtained from three-point QCD sum rule estimates of Refs. [441 145] . 

We should note that in a recent interesting paper [TS] a renormalon analysis of pBj 
and GLS sum rules was also carried out. The aim in that case was to perform a sophisticated 
matching of perturbative and non-perturbative effects in the renormalon subtracted approach. 
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